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ABSTRACT 

We present a mid-infrared high spectral resolution spectrum of CRL618 in the frequency ranges 
778 - 784 and 1227 - 1249 cmT^ (8.01 - 8.15 and 12.75 - 12.85 Atm) taken with the Texas Echclon- 
cross-Echelle Spectrograph (TEXES) and the Infrared Telescope Facility (IRTF). We have identified 
more than 170 ro-vibrational lines arising from C2H2, HCN, C4H2, and C6H2. We have found no 
unmistakable trace of C8H2. The line profiles display a complex structure suggesting the presence 
of polyacetylenes in several components of the circumstellar envelope (CSE). We derive total column 
densities of 2.5 x 10^^ 3.1 x 10^^ 2.1 x 10^'^, 9.3 x 10^^ cm-^ and < 5 x lO^^ cm^^ for HCN, 
C2H2, C4H2, C6H2, and C8H2, respectively. The observations indicate that both the rotational and 
vibrational temperatures in the innermost CSE depend on the molecule, varying from 100 to 350 K for 
the rotational temperatures and 100 to 500 K for the vibrational temperatures. Our results support 
a chemistry in the innermost CSE based on radical-neutral reactions triggered by the intense UV 
radiation field. 

Subject headings: line: identification — line: profiles — surveys — stars: AGB and post-AGB — 
stars: carbon — stars: individual (CRL618) 



1. INTRODUCTION 

The protoplanetary nebula stage (PPN) is one 
of the shortest of a sun-like star's evolution (e.g., 
Ilben fc Renzinil [1983D . Throughout this phase, roughly 
half of the stellar photosphere is ejected, shocking the gas 
of the circumstellar en velope ( CSE) formed in the AGB 
stage (AGB-CSE: Kwok et al.l[r978i) . and unveiling the 
outermost layers of the nucleus. The stellar UV radiation 
field is extremely intense in this phase, photodissociat- 
ing the innermost circumstellar gas and triggering a par- 
ticularly ri ch photochemistry (Woo ds et al.|[2003, 120031 : 
iRedman e t al. 2003; Cernicharo 200i, hereafter C04). 

CRL618 (Westbrook Nebula , AFGL 618) is a very 
young PPN (age ~ 200 yr; iKwok fc BignelllfTgsl . lo- 
cated at a distance of ^ 0.9 — 1.8 kpc f rom the Sun 
jSchmidt fc Cohen' IgSl"; 'Goodrich' 1991 iKnapp et al] 
[1993; Sanchez-Contreras fc Sahai 2004a). It contains a 
BO central star embedded in a dusty ultracompact HII 
region surrounded by a torus and a low-velocity ex- 
panding envelope with an external radius > 20", a to- 
tal mass ~ 1 M0, and an expansion velocity Wexp — 
18.0 - 21.5 km s~| (IKnapp fc Morrislll985l: iFuente et all 
119981: iPardo et al.l |200j). In additi on, it displays gas 
with velocities as high as 200 km s~^ (iBur ton fc Geballd 
[l98l iCernicharo et al.|[T989l: iGammieet al. 1989). This 
high velocity gas (HVG) is the molecular counterpart 
of the brig ht optical jets oriented in the E-W direction 
(|Trammerf '2000) which impact the AGB-CSE and pro- 
duce the well-known opti cal lobes. 

Since its discovery by iWestbrook et al.l (|1975D many 

^ Visiting Astronomer at the Infrared Telescope Facility, which 
is operated by the University of Hawaii under contract from the 
National Aeronautics and Space Administration. 



molecular species have been detected towards this PPN 
(ICer nicharo et al. 1989, 2001a, b; Hcrpin fc Cer n ichar g 
2OOOI : iRemiian et al.l 120051 : ITruong-Bach et al.l Il996[ )7 
of them for the first time in a C-rich 
(e.g., formaldehyde, polyacety lenes C4H2 and 
benz ene, H2O and OH: ICernicharo et al.l 
[2001allbl : IHerpin fc Cernichard 12000^ The 
developed in the AGB phase has been stud- 
in great detail by several authors (see, e.g^ 
2001a, b; Sanchez-Contreras fc SahaJ 
,2004b: Pardo et al 
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2007b, and references therein). 

~!Woods et al. (200l, hereafter W03) and C04 have 
modelled the chemistry of CRL618 suggesting that in 
the innermost envelope the UV photons photodissoci- 
ate C2H2 producing the radical C2H, which can re- 
act with C2H2 to form C4H2 or with H2 reforming 
C2H2. Additionally, C2H can react with C4H2 forming 
C6H2. These processes lead to a rapid C2H2 polymer- 
ization in long carbon chains and clusters. The abun- 
dance ratio between consecutive polyacetylenes (C2nH2, 



1,2,3, 



in CRL618 has been estimated as a factor 



-2-3 (ICernicharo et al.ll2001al hereafter COla). Poly- 
acetylenes are symmetric molecules without a permanent 
dipole moment and, hence, detectable only through their 
ro-vibrational spectra. The strongest bands of their spec- 
tra are expected to fall in the mid-infrared range due 
to the physical conditions prevailing in the innermost 
CSE where the polyacetylenes are built up. However, 
the large telluric opacity in the infrared has largely pre- 
vented the exploitation of this frequency range. This 
issue has been overcome with the launching of the In- 
frared Space Observatory (ISO) and the Spitzer Space 
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Telescope, and the developing of instruments with high 
spectral resolving power such as th e Texas Echelon - 
cross-Echelle Spectrograph (TEXES: iLacv et alj [20021 ) . 
The analysis of these kinds of observations will allow 
us to improve our knowledge about the polymerization 
processes and the formation of complex molecules such 
as long carbon chains, PAHs, and fullerenes (Ceo and 
C70), which have been recently observe d towards circum- 
stellar and interstellar environments (fCami 'eFall [20101: 
IGarcia-Hernandez et all 120101: ISellgren et all 120101) and 
whose ions could be the carriers of the Diffus e Interstel- 
lar Bands (DIBs: lfbmg fc EhrenfremldlfTOOl . 

In this paper, we present a high resolution mid-infrared 
spectrum of CRL618. We have identified ~ 170 ro- 
vibrational lines of bands + vq + + i/g — vg, 
+ + 2z/9 — 2h'g, vq + + v-j — v-j of C4II2, and 
v\\ of C6II2. In addition, we have observed several 
lines of C2H2 and HCN. These lines have been analyzed 
by using a modified version of the mod el of AGE en- 
velopes developed bv lFonfria et alJ ()2008l hereafter F08). 
The observations and the spectroscopic laboratory data 
of HCN and polyacetylenes required to analyze the ob- 
servations are presented in §12] and [31 respectively. A 
description of the spectrum and the model, the adopted 
fitting strategy, and a discussion about the uncertain- 
ties of the parameters can be found in 12) The results 
derived from our fits are presented and discussed in fj5] 
and, finally, summarized in ^Sl 

2. OBSERVATIONS 

We observed CR L618 using the TEXES spectrograph 
(jLacv et all I2002D on the IRTF 3 m telescope during 
the first half of the nights of 12 Dec 2002 and 8 Jan 
2004. TEXES was used in high spectral resolution mode 
which gives R ~ 10^. Six separate settings were used to 
cover from 778 to 784 cm^^ (2002/12/12) and 1205 to 
1250 cm-i (2004/01/08; see Fig. [Tl). For each setting, the 
Beck lin-Neugebauer Object (^BN: iBecklin fc Neugebaueil 
[19671) was observed as a telluric reference before changing 
the instrumental setup. To date, TEXES has detected no 
spectral features in BN at these wavelengths. The spec- 
trograph slit was roughly 1'.'4 wide for all observations. 
The spectral coverage and slit length varied depending 
on the setting. The slit was always long enough that we 
nodded along the slit every 16 seconds to remove back- 
ground emission. 

The data we re processed usin g the TEXES data reduc- 
tion pipeline (jLacv et al.l 120021) . The pipeline operates 
on the raw data files and produces optimally extracted, 
ID spectra with a frequency scale set by user identifi- 
cation of atmospheric lines. Using the difference of an 
ambient temperature blackbody and the night sky emis- 
sion, it is possible to correct partially for telluric spec- 
tral features. Additional correction comes from compar- 
ison with the telluric standard. The baseline has been 
removed by using fourth-order polynomials. The range 
1205 — 1227 cm^^ shows a ripple whose period is roughly 
twice the observed linewidths and that would severely 
compromise analysis of any line in this range. Hence, 
this range has not been analyzed. 

A realistic modeling of molecular lines requires an es- 
timation of the dust emission throughout the envelope. 
It provides us with an estimation of the dust opacity and 
temperature in different shells, allowing us to determine 



how it affects the molecular emission. Since the infrared 
continuum of CR L618 has remained roughly constant for 
sever al decades (jKleinmann et al.l Il978t iPottasch et al.l 
[1981 COla), we have supplemented the TEXES/IRTF 
data with a low resolution SWS/ISO spectrum (COla) 
covering the range between 2.4 and 45 /im. The proper- 
ties of the dusty CSE derived from modeling this contin- 
uum will be subsequently used to calculate the synthetic 
profiles of the molecular lines (see !|4]) . 

3. SPECTROSCOPIC DATA 

All the molecular species considered in this work 
(HCN, C2H2, C4H2, C6H2, and maybe C8H2) are lin- 
ear. The frequencies and intensities of the lines of C2H2 
and HCN have been extensively studied in the laboratory 
but those of C4H2, C6H2, and C8H2 are still not accu- 
rately measured. The frequencies of the ro-vibrational 
transitions of the bands observed in our spectrum are 
well known in the case of C4H2 and C6H2. However, 
the ro-vibrational s pectrum of C8H2 remains unknown 
(|Shindo et al.l[200l . 

C4H2 is a 6-atom linear molecule with nine vibra- 
tional modes. Four of them (y^ — i^g) are doubly- 
degenerated bending modes with energies of 625.4986, 
482.7078, 627.8958 and 219.9778 cm~\ The ro- 
tational constan t in the ground vibra tional state is 
0.1464102 cm-i (jCuelachvili et al.l [198^. The sp ectro- 
scopic constants have been taken from iGu elaclivili et al.l 
(11984"):'McNaughton et al.' ('1 992D : IArie et al.l(|199^ . and 
[Ma tsumur a et al. (1982, 1981). The measured band in- 
tensity is ~ 171 cm- 2 atm-i (jKhlifi et al.lll995D . imply- 
ing that the derived dipole moments of the ro-vibrational 
transitions of band i/g + I's are ~ 0.096 D. 

C6H2 is a 8-atom molecule with 13 vibrational modes; 
the last six (z/g — 1/13) are doubly-degenerated bend- 
ing modes. Their energies are 622.38, 491.00, 258.00, 
621.34, 443.50 and 105.04 cm"!. The rotational con- 
stant in the ground vibrational state is 0.044171 cm~^. 
All these data and the rest of the spectroscopic con- 
stants have been taken f rom [M atsum ura ct al. (1993); 
IMcNaughton et al.l (|1991l ). and IHaas et al.l ([1994). The 
dipole moments of the ro-vibrational transitions of band 
t's + I'll arc ~ 0.1 2 D, derived from a band intensity 
~ 210 cm-2 atm-i (jShindo et al.|[2003l ). 

The central frequency of the C8H2 vio + vu band is 
not accurately known but it is expected to fall inside the 
observed range at ~ 1230 cm"^. It h as a small ground ro - 
tational constant, ~ 2 x 10"^ cm'^ (jShindo et al.ll200lD . 
which will produce a band formed by a series of lines 
partially resolved with the resolution power provided by 
TEXES, and a FWHM ~ 5 cm'^ at the temperatures 
prevailing in the innermost CSE of CRL618. This struc- 
ture could be easily identified if the abundance of this 
species is large enough to produce absorption above 10% 
of the continuum, in view of the large density of lines 
from the other species. The estimated dipole moment of 
this band is ~ 0.14 D, which has b een derived from th e 
band intensity ~ 256 cm'^ atm'^ CShind o et 311 12001'). 

For these three molecules (C4H2, C6H2, and C8H2), 
only the intensities of the whole bands are available in 
the literature. It is worth noting that the observed bands 
are blended with several hot bands (e.g., for C4H2 the ob- 
served band is Uq + and the associated hot bands are 
V6 + 1^8 + - ^9, + + 2i^9 - 2vq,. . . ) involving vi- 
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Fig. 1. — Observed and synthetic spectra towards CRL618 (black and red), (a) Spectra in the range 778 — 784 cm ^. The detected 
hnes are produced by C2H2 and HCN (green and blue). The asterisks (yellow) indicate several absorptions that could be assigned to the 



Hi^CCH lines i/5Re(21), Re(22), and Re(23) (* 



and 



respectively), (fe and c) Spectra ranging from 1227 to 1249 cm contain 



C4H2 lye + us and C6H2 lyg + '^ii bands (green and blue). See the text for the modeling details of these lines. [See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 2.— Observed continuum of CRL618 (SWS/ISO; black) and 
fit between 6 and 45 /im (red). A broad feature can be seen at 
~ 28 fim, which can be assigned to the solid state band produced 
by the unknown material usually assumed to be MgS. The fit has 
been calculated by assuming dust grains composed of amorphous 
carbon (see text for details). 

brational states significantly populated at room temper- 
ature. Therefore, we have estimated the dipole moment 
of the ro- vibrational transitions of the fundamental band 
of each molecule from the intensity measurements quoted 
in the literature and we have assumed that the hot bands 
{vq + 1/^ + 1^9 — 1^9, 1^6 + 1^8 + 21^9 — ^vg,. . . ) have the same 
dipole moment as the vq + one. The consequences of 
this approximation will be discussed in tj4.1l 

The frequencies and dipole moments of C2H2 and 
HCN have been tak en from the HITRAN Database 
(jRothman et al.ll2005f ). 

Several bands from ^'^C-bearing isotopologues of HCN 
and the polyacetylenes a re expected to fall in the ob- 
served frequency ranges. iPardo fc Cernicharo I ()2007al ) 
suggested a ^'^C/^'^C ratio ~ 15 in the torus, mean- 
ing that [HCN]/[H^3CN] ~ 15 and [CzHaJ/p^CCH] ~ 
[C4H2]/[Hi^CCCCH] ~ [C4H2]/[HC"CCCH] ~ 7.5, 
where H^^CN, H^^CCH, H^^CCCCH, and HC^^CCCH 
are the most abundant isotopologues after the main ones. 
The HCN lines present in the observed spectrum are 
too weak to expect a detectable absorption from H^'^CN 
(Fig.[T]a). Concerning acetylene, although H^'^CCH lines 
might be assigned to several very weak features (Fig. [Da), 
their opacities are too low to affect the much stronger 
C2H2 and HCN lines or to be used to derive any reliable 
information about H^^CCH. There is almost no informa- 
tion about the isotopologues of C4H2 in the literature. 
The center of the bands have been estimated from ab 
initio calculations to be shifted betw een 5 and 10 cm~^ 
from that of the main isotopologue ()McNaughton et aLl 
[1992), but these values could be affected by large un- 
certainties. All the significant features in the frequency 
range 1228 — 1249 cm~^ of the observed spectrum can be 
explained by considering just C4H2 and C6H2 (Fig. [3) 
and c). Therefore, either the bands of the isotopologues 
actually fall out of the observed frequency range or their 
abundances are insufficient to produce detectable lines. 

4. DESCRIPTION OF THE SPECTRA AND MODELING 

The broad continuum of CRL618 taken with SWS/ISO 
displays maximum emission at ~ 40 /im, typical of a 
cold envelope (see Fig. [2]). The spectrum shows solid 
state bands, although much weaker than those displayed 
in the continuum of the C-rich AGB star IRC-l-10216, 
where the bands due to solid SiC at 11 /im and that at 
27 jjLVD. (whose carr ier is still unkno wn but usually as- 
signed to MgS; e.g..lHonv et al.ll2002l ) are clearly present 
(jCernicharo et al.lll999() . Despite the complex geometry 



found in the CSE of CRL618 (see fJD), the broad con- 
tinuum emission from the dust grains smooths any ge- 
ometrical effect in the spectrum. Therefore, the dusty 
AGB-CSE can be assumed to be spherically symmetric 
in a first approximation. We have used the model of the 
CSE developed by F08 and assumed t hat the dust grains 
are c omposed of amorphous carbon (jRouleau fc Mart"inl 
|1991|) to fit the continuum between 6 and 45 /xm. To 
date, this model cannot deal with radiative scattering 
by dust grains. Hence, we limit our fitting to wave- 
lengths above 6 /im. The central source has been as- 
sumed to have an angular size Ur ~ 0'.'135 (jPardo et al.l 
[200ll2007bn . which was determined from a high-quality 
modeling of a large set of HC3N lines in the millimet- 
ric range. We have adopted a distance to the source of 
900 pc, derived from the luminosity-distance relationship 
for PPNe and estimations of the proper motion of dense 
inhomogeneities in the high velocity lo bes (Goodrichl 
119911: isinchez-Contreras fc Sahail l2004a[) . These data 
imply a linear radius for the continuum source Rc — 
1.8 X 10^^ cm. The CSE is assumed to have a radius 
of 600 ac, since the model displays no significant contri- 
bution to the fiux from gas and dust located at larger 
radii. We have chosen a dust temperature following a 
r"^"* law with different values for the exponent in the 
torus and the AGB-CSE. The exponent 7^ has been as- 
sumed to be ~ 0.39 over the latter regfo n as in the case 
of IRC-l-10216 (F08; [SocEiEaD [1983). The synthetic 
spectrum has been fit to the data by eye, the best fitting 
strategy when features such as solid state bands, which 
cannot usually be well reproduced with current labora- 
tory or theoretical data, are present. The parameters 
derived following this procedure are shown in Table [TJ 
The largest differences between the synthetic and the ob- 
served spectra come from the C2H2 and HCN bands at 
13 /im and the solid state features between 20 — 40 /im 
which are not considered in the fit (Fig. [5]) . A discussion 
about the sensitivity of the synthetic spectrum to varia- 
tions in the parameters presented above can be found in 

We estimate the rms noise of the TEXES/IRTF spec- 
trum to be cr ~ 1%. With this noise level, we assume 
that a line is detected if its absorption and/or emis- 
sion is larger than 3tT. The part of the spectrum rang- 
ing from 778 to 784 cm~^ shows three C2H2 lines from 
band 1^5 and two HCN lines from band V2- However, 
there are some features in the spectrum that do not 
fulfil the 3ct condition but might be assigned to tran- 
sitions of the C2H2 bands 1^4 -I- 2^5 — i^4 and 2i^5 — 1^5. 
The strongest line in this range is C2H2 i^5Re(21). The 
range between 1228 and 1249 cm~^ shows C4H2 lines 
from the fundamental band vq + vg, and the hot bands 
v& + V?, + 1^9 — vq, vq + Us + 2i^9 — 2i^9, and (maybe) 

+ i/g + 1^7 — V7 (ro- vibrational co nstants for these bands 
are from iMcNaughton et al.iri992| ). The v^ + vu band of 
C6H2 appears as weak lines often overlapped with those 
of the C4H2 + v^Pe branch (Figs. [Hand |3]). 

Both C2H2 1^5 and C4H2 I'e + lines are broad 
(FWHM ~ 35 and 20 km s^^ with a spectral resolu- 
tion of ~ 4 km s^^), in contrast with those of the hot 
bands (c± 10 km s~^; Fig. [3]). The profiles of the lines 
of these fundamental bands exhibit a main absorption 
surrounded by several features at lower and higher fre- 
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TABLE 1 

Derived Parameters 



Parameter 


Units 


Value 


Errors 


IN I Xl.V_ylN 1 




(0 ^ 4- Q\ \^ ^ n^'^ 


I"; 


NCCoHo) 


_2 

cm 


(3.1 ± 0.7) X 10^'^ 


(2) 






?2 1 + ■^"1 X in^'^ 






_2 

cm 






Vexp(Kc) 


km s 


1 nn+0-40 
1-00-0.22 


(2) 


Vexp{Rdl) 


km 


8-00t°-J° 


(2) 


VexpiRd2) 


km s^^ 


17.0 ±0.8 




R-dS 


R 


11.0 ±1.0 




T„t{HCN) 


K 


350 ± 50 


(6) 


Trot(C2H2) 


K 


200 ± 20 




Trot(C4H2) 


K 


100 ± 20 


8! 


^rot(C6H2) 


K 




(5) 




K 


400 ± 50 


(1) 


TdiRc) 


K 


110 ±6 




TdiRds) 


K 


98 ±4 


\i 


Td{8 fim) 




2.3 ±0.3 




^torus 


% 




(1) 


Cage 


% 


35 ± 5 


(1) 



Note. — N(X): total column density of species X; Vcxp(t): 
gas expansion velocity at radius r; Rds '- radius of the outer bound- 
ary of the torus; T^otCAT): rotational temperature of species X at 
the inner boundary of the torus; Tc'. black-body temperature of 
the central source; Td{r): dust temperature at radius r; Td{8 fj,m): 
dust optical depth at 8 ^m; ^torus ^-i^d CagB ■ fraction of the total 
dust optical depth coming from the torus and the AGB-CSE. The 
uncertainties have been calculated by varying the parameters re- 
lated to (1) continuum, (2) C2H2 !^5Re{21), (3) C4H2 U(i+usRe{7), 
(4) C4H2 U6 + i^8Re{22), (5) C6H2 us + i^iiRe(34), and (6) HCN 
!^2Re{22). N(C2H2) and N(HCN) have been estimated from the 
CSE and the HVG at once. See J5] for separate estimations in 
each structure. In the calculations, we have assumed the distance 
t o CR L618 and the size of the central source as fixed. See the text 
(" ^14.11 1 for a deeper explanation about the strategy followed for the 
estimation of the uncertainties. 



quencies. Some of these features are associated with 
hot bands and with CrH q hues from t he fu ndamental 
band i^s + i^ii (Fig. E]). IPardo et all (|2004 hereafter 
P04) suggest a constant turbulence velocity in the torus 
~ 3.5 km s^^ and an expansion velocity < 18 km s^^, 
allowing many of the unassigned features within the line 
profiles to be attributed to C4H2 and C2H2 expanding 
in this region. However, the origin of the remaining 
C2H2 line features is less certain since they require ex- 
pansion velocities significantly larger than that observed 
in the torus (see Fig. [3]). We estimate from our model 
that the emission from the central source at 8 fim is 
at least five orders of magnitude larger than that com- 
ing from regions with an impact parameter b > ac- 
Hence, these unexplained features could arise from gas 
in the HVG located in front of the central source (P04; 
iSanchez-Contreras &: Sa hai 2004a). 

The model adopted to interpret the TEXES/IRTF 
observations is a modified version of that devel- 
oped by F08. Although the innermost CSE displays 
roughly axial symmetry, as several interferometric obser- 
vations esta blished (Sanch ez-Contreras fc Sahail l2004al : 
ISanchez-Contreras et al.ll2004bf) . we have assumed spher- 
ically symmetric abundance distributions for the consid- 
ered species as a first approximation. This is a reason- 
able approach since the shells near the central source 
should not be significantly affected by the outflows. The 
CSE has been divided in four regions {Zi, . . . , Z4) re- 
lated to the torus and the AGB-CSE (see Fig. g]). We 
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Fig. 3. — Comparison between observed and synthetic C2H2 and 
C4H2 lines. The strongest lines are (a) C4H2 + i^sC'^'u )Re(19), 
(b) C4H2 V6 + i^8(o-J)Pe(21), and (c) C2H2 i/5(7r„)Re(21). These 
spectra show additional C4H2 and C6H2 lines (panel a: Li and 
L2 are C4H2 + i^s(fu )R'e(20) and Re (18); panel 6: C4H2 in 
blue and C6H2 in green). The velocity axis is centered in the 
LSR systemic velocity (~ -21.3 km s"!; IPardo et al.l[200D . The 
maximum absorptions in the strongest lines arise from the CSE, 
displaying velocities of ~ —3.5, — 9.0, and —16.0 km s~^. The 
C2H2 line (panel c) shows the HVG contribution. These features 
are not present in the C4H2 lines implying that they arise just 
from the CSE. [See the electronic edition of the Journal for a color 
version of this figure.] 



have assumed a gas density profile oc 



cxp 



(constant 



mass-loss rate), excitation temperature profiles (vibra- 
tional and rotational) oc r~", with a depending on the 
region, and initially a co nstant turbulence velocity of 
3.5 km s^i in the torus (jPardo et al.ll2004D . We have 
added three gas condensations in the inner boundary of 
region Zi with different physical and chemical conditions 
with respect to the innermost CSE in order to model an 
inhomogeneous HVG. The synthetic molecular lines have 
also been fit to the observed ones by eye. This strategy 
allow us to fit lines partially overlapped with unknown 
features or affected by irregular baselines impossible to 
be removed (see Fig. [Ij)). The derived parameters are 
shown in Table [T] 

4.1. Uncertainties 

We estimated our uncertainties as in F08. Specifically, 
we explore how the variation of a given parameter affects 
either the continuum or a molecular line, depending on 
the nature of the parameter. When using a molecular 
line, we choose the one most sensitive to the variation of 
the parameter. We consider a synthetic continuum/line 
as a good fit to the observed one when all its points devi- 
ate from those of the best fit in less than 10% the maxi- 
mum absorption when modeling a molecular line and the 
maximum emission when modeling the continuum. Oth- 
erwise, the synthetic spectrum is considered as a bad fit. 
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Fig. 4.— Model of the inner envelope of CRL618. The CSE is 
composed of four shells: Z\, and Z-^ inside of the torus, and 
^4 corresponding to the AGB-CSE. We have adopted an angu- 
lar radius for the central source clc — 0''135 {P04). The values 
of the angular radii a^i and 0^2 are assumed to be 1.5ac and 
2.0ac, respectively, to allow the model to deal with unexpected 
sharp variations of the physical conditions. The radius a^^ has 
been considered as a free parameter. We have assumed a constant 
turbulent velocity tit^j. ~ 3.5 km s~^ in the torus (P04). The high 
velocity gas (HVG) has been modeled as gas condensations (see 
text). The axis on the left shows the angular offset with respect 
to the position of the source. On the right, we have included the 
point spread funcion (PSF) of the telescope to allow comparisons 
with the model. 

The minimum and maximum values allowed for a given 
parameter are the smaller and larger of its values which 
produce good fits when varying, at the same time, all 
the other parameters. Therefore, the lower and upper 
uncertainties are the differences between the minimum 
and maximum values and that derived from the best fit. 

In this work, we have fixed two parameters that cannot 
be estimated from our observations at the same time as 
the rest: the distance to CRL618, D, and the angular 
radius of the central source, ac- To date, these parame- 
ters are not well known, specially D, which could be two 
times larger (see SjlJ. 

A variation in any of these parameters affects the lin- 
ear radius of the central source. An increase in the acD 
product would raise the optical depth {ti, oc acD) dimin- 
ishing the escape probability of the infrared photons that 
come from the the inner shells of the torus. This effect 
could be counteracted by decreasing the column density 
of the considered molecular species and the dust up to a 
factor ~ 2 — 3 in our case. 

In addition, a variation in ac would have an impact on 
the temperature profiles as well. The most affected pa- 
rameters would be the exponents while the temperatures 
at the inner boundary of the torus would require minor 
corrections to recover the best fit. 

The approximations adopted to estimate the dipole 
moments of the ro- vibrational transitions of C4H2, C6H2, 
and C8II2 also introduce additional uncertainties in sev- 
eral parameters. The dipole moment of a ro-vibrational 
transition behaves as a second degree poly nomial on 
J(J -|- 1) (see, e.g., iJacquemart et all l200l( i. while we 
have assumed that all the ro-vibrational transitions in a 
band have the same dipole moment. This different de- 
pendence on J alters the derived column densities and 
excitation temperatures from the real ones. Adopting 
a maximum deviation of 10% for the (constant) dipole 
moment squared used in our model with respect to the 



real one (as occurs for C2H2; IJacquemart et all 120011) . 
the derived column densities and excitation temperatures 
would be affected by an additional error up to 30%. 

Similar variations in the uncertainties of the vibra- 
tional temperatures could occur due to having assumed 
no differences between the vibrational dipole moments of 
the fundamental band and the associated hot bands. 

5. RESULTS AND DISCUSSION 

The fits indicate that the deepest absorptions of 
all lines (Fig. [3]) have velocities ~ —3.5, —9.0, and 
— 16.0 km s~^, and come from regions Zi, Z2, and 
respectively. This is compatible with an expansion ve- 
locity profile in which Wcxp(^c) — 1-0, Wcxp(^di) = 8.0, 
Wexp(-Rd2) = 17.0, and Woxp(?' > Rds) = 18.5 km s"^ The 
line C2H2 j^5Re(21) also displays several features with ve- 
locities ~ 10.0, 2.3, and —21.8 km s~^, not present in the 
C4II2 lines (Fig. [3]). They do not arise from the CSE and 
are not hot bands. This fact supports the existence of 
C2H2 in the HVG and suggests that the abundance of 
C4II2 is negligible in it. 

The discrete absorptions found in the C2II2 lines 
(Fig. [3]) produced by the HVG could be the fingerprints 
of dense clumps whose outermost layers protect the inner 
molecular gas from being dissociated by the UV radiation 
field impinging on them and allow a sig nificantly slowe r 
chemical evolution in their cores (Re dman et al.|[2003() . 
In these circumstances, these clumps must be formed due 
to («) an irregular mass- loss in the latest phas es of the 
stellar atmosphere ejection (jPvson et al.lll989| ). or {ii) 
density fiu ctuat i ons or instabilities in th e ioniza tion front 
JCaDrfotti 1973HGarcfa-Segura k Fran co 199(1: iWilham^ 
[l999i) . As the clumps must be in the HVG and in front of 
the central source, we estimate that their characteristic 
size should be < 1 i?c. 

We derive a vibrational temperature, Tvib, — 350 K 
for C2H2 in the innermost CSE and in the HVG. In the 
CSE, it decreases to ~ 200 K between Rc and Rdi reach- 
ing ~ 50 - 100 K at Rd3. For C4H2, Tvib ^ 500 K in the 
innermost CSE followed by a fast decrease to ~ 400 K at 
Rdi and reaching ~ 200 — 250 K at the end of region Z3. 
It is not possible to derive a reliable Tvib for C6H2 since 
its lines are very weak but the fits support values similar 
to that of C4H2. The HCN lines do not allow an accurate 
Tvib determination in the torus or in the HVG because 
of its low abundance. However, the lack of absorption 
from the hot bands and emission from the fundamen- 
tal band implies Tvib ^ 100 K. The difference in Tvib 
between HCN, C2H2, and C4H2 is likely produced by 
different pumping mechanisms between these species. In 
particular we note that some of the hot bands arise from 
metastable vibrational states (i.e., infrared- inactive) . 

The synthetic spectra show that the C2H2 and 
C4H2 hot bands require lower turbulence velocities (< 
2 km s""'^) than fundamental bands (~ 3.5 km s~^). It 
suggests that the turbulence velocity would be lower in 
regions Zi and Z2, where most of the hot bands arise, 
than in the rest of the torus. We would need additional 
observations with a higher signal-to-noise ratio and a 
more realistic model to explain this result. 

Concerning the rotational temperatures, Tjot, we de- 
rive Trot — 200 K for C2H2. This value has little error 
because the observed lines are very sensitive to variations 
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in Trot- For C4H2, we derive Tj-ot — 100 K also accurately 
due to the large number of observed lines. In the case of 
HCN, Trot — 350 K but this time the uncertainties are 
large since we have only observed two weak lines. The 
diversity of rotational temperatures indicate that these 
species are out of LTE, at least in regions Z\ and Z^^ 
which dominate the molecular absorption. 

The results support a 3:1 ortho-para ratio for the poly- 
acetylenes, i.e., the ratio expected under LTE conditions. 
This fact implies that the formation of polyacetylenes 
occurs in warm regions as suggested in previous works 
(e.g., COla; C04), where the ortho-para ratio of the 
new molecules is given by the spin degeneracy. Taking 
into account that the rotational contants of the poly- 
acetylenes are small compared to the kinetic tempera- 
tures prevailing in the emitting regions, we do not expect 
ortho-para conversions even though they are not strictly 
forbid den, contrarily to the case of H2 at low tempera- 
tures (|Herzberelll963D . Moreover, our result also points 
out that other processes capable of producing ortho-para 
conversions such as (?) proton exchange with atomic hy- 
drogen, (ii) collisions with paramagnetic molecules or 
ions such as H+ , , or H;^ , and {iii) spin inversion pro- 
cesses on the dust grain surfaces (lFarkaslll935l: iHerzberd 
[ToellMassie fc Huntenll 1981 Burton et al.lll992n . are in- 
efficient in the innermost CSE. 

The C2H2 column density in the CSE is ~ 2.0 x 
10^^ cm~^ coming 45% from region Zi, 30% from Z2, 
and 25% from Z3. The contribution from the AGB-CSE 
(region Z4) is negligible. The averaged column density 
in the HVG is ~ 1.1 x 10^'' cm^^. Hence, the total col- 
umn density is ~ 3.1 x lO^'' cm~^, in good agreement 
with 2 X 10^^ cm^^, the value derived by COla from 
their low spectral resolution SWS/ISO observations (see 
Table [T] for an error estimation). In addition, the col- 
umn density that we propose for the CSE agrees with re- 
sul ts of the time -dependent chemical models developed 
bv iWoods et al.l ([2002, 2003), which bracket it in the 
range ~ 5 x 10^^ — 5 x 10^® cm^^, depending on the 
evolution degree of the envelope. It is possible to esti- 
mate the C2H2 abundance ratio between the clumps ex- 
panding in the HVG and the gas in the innermos t CSE 
and to compare this result with those of Redman et al.l 
(|2003f) . Assuming a characteristic length in the line of 
sight < 0.5 i?c, the averaged gas density in the clumps 
is > 100 cm"'^. Hence, as the C2H2 density in the inner- 
most CSE is 300 cm~^, th e ratio is > 0.3. Following 
the results bv iRedman et al.l (|2003[ ). this lower limit in 
the ratio would be achieved after ~ 1100 yr of evolution. 
However, the age of CRL618 as a PP N has been esti- 
mate d in several hundreds of years (e.g.. lKwok fc BignelH 
[l98l . 

The C4H2 column density in the CSE is ~ 2.1 x 
10^^ cm-2 (35% from region Zi, 60% from Z2, and 5% 
from Z3), again with insignificant contribution from the 
AGB-CSE (Table As we do not see C4H2 in the 
HVG, the total column density is that derived for the 
CSE which is a factor ~ 2 larger than 1.2 x 10^^ cm~^, 
suggested by COla. The difference between both values 
is probably due to the large PSF of ISO compared to that 
of IRTF (~ 3" and 0'.'7 at ~ 8 ^m, respectively). In this 
case, our result is very similar to 2.3 x 10^^ cm~^, the 
maximum abundance for C4H2 according to the model of 



W03 for the evolution of the CSE. As the abundance of 
C4H2 in the HVG is much lower than in the innermost 
CSE, we estimate the C4H2 abundance ratio between 
the clumps and the innermost CSE in a factor < 0.1. 
This upper limit supports the results by IRedman et al.l 
([2003). which suggest a ratio ~ 0.01 — 0.3 for clumps 
aged between 100 and 1000 yr. This period of time and 
that derived from the C2H2 results suggest that a faster 
chemistry could be needed to better reproduce our re- 
sults. However, the discrepancies between their results 
and ours could have to do with differences in the gas 
density, directly related to the adopted distance to the 
star, or different simplifications regarding the detailed 
physical structure of PPN (see more complex model s in 
iPardo et al]|2004l: iSanchez-Contreras fc Sahaill2004al) . 

For C6H2, the column density in the CSE is ~ 9.3 x 
lO'^^ cm~^ (15% from region Z\ and 85% from Z2; see 
Table [1]) . There is no detectable contribution to the ob- 
served C6H2 lines from region Z3 and from the HVG. 
Hence, the total column density is in good agreement 
with the value of 6 x 10^^ cm~^, derived by COla. The 
maximum column density proposed by W03 (~ 7.7 x 
10^^ cm~^) is similar to our result but somewhat lower. 
It might indicate that the model developed by these au- 
thors needs to be slightly improved as far as the inter- 
action between polyacetylenes and dust grains or uncon- 
sidered chemical reactions involving polyacetylenes are 
concerned. 

Finally, the column density of HCN in the CSE is 
~ 2.0 X 10^^ cm~^ and the averaged column density in 
the HVG is ~ 5.0 x 10^^ cm^^^ g^^jj results are affected 
by a large error due to the weakness of the observed lines. 
Hence, we can conclude that the total column density is 
~ 2.5 x 10^^ cm~^ (Table [T]), in good agreement with 
1.5 X 10^^ cm~^ proposed by COla but sma ller than the 
value of ~ 4 — 7 x 10^^ cm~^, derived bv IPardo et al.l 
(|2004 I2005D . The discrepancy existing with the lat- 
ter result can be attributed to HCN molecules known 
to exist in the cold region Z4 (AGB-CSE), unobserved 
through infrared observations. As in the case of C4H2, 
the maximum column density expected by W03 is sim- 
ilar but larger than that derived by us, suggesting that 
their model works also fine for HCN. 

The abundance of C8H2 estimated from the results by 
COla and that proposed by W03 and C04 suggests that 
the C8H2 v\Q -\- 1^14 band should be observed in our spec- 
trum. Nevertheless, we have not found any pattern of 
lines that could be assigned to this species. The tel- 
luric contributions to the raw spectra of CRL618 were 
removed by dividing the latter by the corresponding spec- 
tra of the BN object, free of intrinsic molecular lines in 
the observed frequency range. Further corrections were 
performed to mend minor instrumental features from the 
spectra. Therefore, we are sure we did not accidentally 
remove any broad feature such as the C8H2 band. The 
upper limit to its total column density, calculated by 
assuming a peak intensity for the v\q -f 1^14 band lower 
than 10% the continuum, and the same excitation tem- 
peratures and coexisting with CgH2, is 5 x 10^^ cm~^. 
This value is lower than the column densities of C2H2, 
C4H2 and C6H2 in the CSE by a factor 4.0, 4.2 and 1.9, 
respectively. 

The abundance ratio [C2H2]/[C4H2] in the torus is 
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~ 0.95, compatible with a region where the abundances 
have reached the steady state (C04). The abundance ra- 
tios in each region {Zi, Z2, and Z3) are 1.2, 0.49, and 3.5, 
respectively. The abundances in Zi seem to be in steady 
state while in Z^ they are still evolving. In region Z2 the 
ratio is very low indicating variations in the abundances, 
likely related to condensation onto dust grains or involve- 
ment in chemical reactions still unconsidered in the cur- 
rent chemical models. The ratio [C4H2]/[C6H2] ~ 2.3, 
in good agreement with the results suggested by C04 for 
evolving abundances. The C6H2 lines are too weak to 
derive reliable ratios for each region in the torus. 

5.1. SMP LMC 11, CRL618, and IRC+10216 

The protoplanetary nebula SMP LMC 11, located in 
the Large Magellanic Cloud, is known to display several 
bands in absorption arising from C2H2, C4H2, and C6H2 
(|Bernard-Salas et al.ll2006r and references therein) . This 
fact suggests that the chemistry in its CSE is very similar 
than in CRL618. Hence, a comparison between these 
two sources could contribute interesting clues to a future 
chemical model of SMP LMC 11. 

The observations carri ed out in the infrared by 
[Bernard Salas et all (|2006[ ) towards SMP LMC 11 show 
a wide absorption in the continuum ranging from 12 
to 17 ^m, formed by the blending of several hundreds 
of lines of C2H2, C4H2, C6H2, and other abundant 
molecules such as CgHg. From this feature can be in- 
ferred that C2H2 is more abundant than in CRL618, 
even as much as in IRC-l-10216 (with a column density 
~ 1.6 X 10^^ cm~^; F08), where a broad absorpt ion at 
~ 13 /J,m was also observed (jCernicharo et allll999[ ). The 
presence of strong hot bands such as va + i^5 — vz, 2^5 — 
(at 13.5 — 14.0 /^m), and v^ + v^ (at ~ 7.5 /im) in absorp- 
tion indicates that the vibrational temperature in the 
innermost CSE should be ~ 500 — 600 K, i.e., between 
those of CRL618 (with weak hot bands; see Fig. [H^a)) 
and IRC+10216 (with strong hot bands; F08). This ex- 
citation would be probably caused by a somewhat strong 
infrared radiation field produced by the warm dust, as in 
the case of IRC-)-10216 (e.g., F08). 

Concerning the rest of the polyacetylenes, the lower 
abundance ratio between C6H2 and C4H2 in SMP LMC 
11 than in CRL618 suggest than the former source is 
chemically less evolved than the lat t er (^y 03,C04), in 
agreement with [Bernard- Salas et al.l ()2006D . To date, 
it has been impossible to clearly detect C8H2 towards 
CRL618. However, the great resemblance in their 
chemistry and the good agreement between observa- 
tions and chemical models suggest that C8H2 could 
be built-up from CfiH2 i n the CSE of both sources. 
iBernard-Salas et al.l (|2006[ ) compare in their Figure 2 two 
low resolution mid-infrared spectra of CRL618 and SMP 
LMC 11. These spectra seem to share a weak unidenti- 
fied feature at ~ 16 ^m, between the bands and vn of 
C4H2 and C6H2, respectively. This feature could be the 
Q branch of the C8H2 band 2^14, se veral times stronge r 
than J^io + i^i4 at room temperature (jShindo et al.ll200ll ). 
It should be observed in the future with higher spectral 
resolution to assess whether it can be assigned to the 
Via band of C8H2 . The recent discovery of Cgo and C70 
(|Cami et al.ll2010D strongly support the growth mecha- 
nism for carbon clusters proposed by C04. In his models, 
carbon clusters where produced from the photodissocia- 



tion of polyacetylenes and, although limited to Cig, the 
predicted abundances for these large carbon species is 
very large. 

6. SUMMARY AND CONCLUSIONS 

In this Paper, we present high-resolution mid-infrared 
observations towards the PPN CRL618. The sampled 
spectral range (778-784 and 1227-1249 cm"^) observed 
with the high resolving power spectrograph TEXES allow 
us to resolve bands v^+vg, and vg+vn of C4H2 and C6H2, 
respectively, in addition to several lines of bands and 
1^2 of C2H2 and HCN. These rich data have enabled the 
modeling of the useful ro- vibrational line profiles of these 
molecular species to estimate their abundances and the 
physical conditions of the gas and the dust throughout 
the inner circumstellar envelope. 

The analysis of the observations has yield the following 
results among others: 

• our results support the chemical model suggested 
by W03 and C04 for the polymerization of C2H2; 

• most of the HCN and C2H2 are in the inner CSE. 
The rest come from several dense clumps located 
in the high velocity gas; 

• C4H2 and C6H2 are formed in the innermost CSE. 
Their abundances seem to be negligible in the 
clumps since the emission from these molecules is 
undetectable in our spectrum; 

• we are not able to detect any trace of the C8H2 
band vio + vu, expected to fall in the observed 
range. This implies an even lower abundance for 
this species compared to previously suggested val- 
ues. An upper limit to its column density has been 
estimated; 

• there exist large differences between the excitation 
temperatures (vibrational and rotational) of HCN, 
C2H2, and C4H2 which indicate that the inner CSE 
is out of LTE. 

In addition, the results of this work demonstrate the 
power of IR observations in the determination of the 
abundances and physical conditions of the gas in com- 
plex structured environments such as the innermost en- 
velopes of the evolved stars. Further improvements in 
the search for C8H2 will be made in the future by ob- 
serving CRL618 and SMP LMC 11 at ~ 16 ^im with the 
Echelon-cross-Echelle Spectrograph (EXES) mounted on 
the Stratospheric Observatory for Infrared Astronomy 
(SOFIA). 
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